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1. Introduction 



The bimetric theory of gravity were introduced in ^ in order to describe the interaction 
of gravity with a massive spin 2— meson. These theories were also investigated recently 
due to their cosmological solutions [Q, |3|, ^, |5|, ^, |^ ^, ^, ^] . Generally, bimetric theories 
contain one massive and one massless spin— 2 field and one scalar ghost mode. As a result 
the bimetric theories of gravity are plagued by the same problems as the general massive 
theories of gravity. 

On the other hand recent formulation of the massive gravity known as the non-linear 



massive gravity seems to be free of ghosts |11, 12 1, for further improvement, see [14, 13 1. 



This theory was further extended in |15| where the theory was formulated with general 
reference metric. The Hamiltonian analysis of given theory was performed in several papers 



m, m, I 



21, 22 1 with the most important results derived in [23, 24] with the 



outcome that this non-linear massive theory possesses one additional constraint and the 
resulting constraint structure is sufficient for the elimination of the ghost degree of freedom. 

Due to the fact that the non-linear massive gravity with general reference metric is 
free of ghost it is tempting to generalize this theory by introducing the kinetic term for 
the reference metric which now becomes dynamical. As a result g^i^ and /^j/ come in 
the symmetric way in the action which means that the non-linear massive gravity was 
generalized to the bimetric theory of gravity. This important step was performed in [|25|. 



Then it was argued in [23[ that the resulting theory is the ghost free formulation of the 
bimetric theory of gravity using the very detailed proof of the absence of ghosts in the 
non-linear massive gravity performed here. However the question of the absence of the 
ghosts in the new formulation of the bimetric theory of gravity was questioned recently in 
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|28] ^. This analysis was based on the careful Hamiltonian analysis of the bimetric gravity 



with redefined shift function |25]. We have argued that due to the fact that the lapse 
function of the metric f^i, is now Lagrange multiplier whose value is determined by the 
dynamics of the theory it is not possible to find another additional constraint so that the 
ghost mode is still present. On the other hand we showed that there are no first class 
constraints corresponding to the generators of the diagonal diffeomorphism and we argued 
that this is a consequence of the redefinition of the shift function performed in |25]. 



The goal of this paper is to extend the analysis presented in to the case of the 
Hamiltonian analysis of the bimetric theory of the gravity with general potential between 
two metric ^^1/ and f^u- We identify four the first class constraints corresponding to the 
diagonal diffeomorphism and determine the Poisson brackets among them in order to show 
that they obey the right form of the algebra. These results immediately solve the issue 
found in since we are now able to identify the four first class constraints that generate 
the diagonal diffeomorphism in case of non-linear bimetric gravity. As the next step in our 
analysis we try to answer the question whether it is possible to find an additional constraint 
in bimetric gravity which could eliminate the ghost mode. In other words we analyze the 
time development of the four remaining constraints. We present some evidence that it is 
very difficult to have such a form of the potential that will lead to the emergence of the 
additional constraint. Then we analyze the square root form of the potential and we argue 
that it is difficult to find an additional constraint which suggests that the ghost mode is 
still presented in given theory which also confirms the observation presented in |28|. 



The structure of this paper is as follows. In the next section (|2|) we introduce the 
bimetric theory of gravity and find its Hamiltonian formulation. Then in section (|3|) we 
calculate the Poisson brackets between constraints and find the constraints that are gen- 
erators of the diagonal diffeomorphism. In section (Q) we analyze the time evolution of 
the second class constraints. Finally in conclusion (|^ we outline our result and suggests 
possible extension of this work. 

2. Hamiltonian Formulation of Bigravity with General Potential 

Let us consider bimetric theory of gravity defined by the action 

S = Mlj d^x^gR{g) + Ml j d^x^/^fRi f)-f, j d\{det g det ffl^ViH^',) , (2.1) 
where we presume that the potential term is the general function of 

= . (2.2) 



When we require that the action (2J) is invariant under following diffeomorphism trans- 
formations 

dx^ dx^ " dx^ dx^ 



""^This issue was also addressed in the recent paper 
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we have to demand that the potential term has to be functions of various powers of i/^ and 
their traces. The goal of this paper is to extend the Hamiltonian analysis of the particular 
bimetric gravity performed in |3C] to the case of the general form of the potential V. 

To begin with we introduce the 3-1-1 decomposition of the four dimensional metric Q^u 

goo = -N"^ + Nig'^Nj , goi = Ni , gij = gij , 



^00 -o^^^ -ij ^ a N'N^ 



(2.4) 



together with the metric /^^ 



foo — — -I- Lif^Lj , foi — Li , fij — fi. 



fOO _ fOi fij fij _ ji T fji 



(2.5) 



To proceed further we use the well known relation ^ 



(2.6) 

where i?^^) and R^^^ are three dimensional scalar curvatures evaluated using the spatial 
metric g^j and fij respectively and where the extrinsic curvatures Kij and Kij are defined 

1 . 1 . . 

Kij = — {dtgij - ViNj - VjNi) , Ki, = — {dtfij - ViLj - V^U) , (2.7) 

and where Vj and Vj are covariant derivatives evaluated using the metric components gij 
and fij respectively. Finally note that Q^^^'- and Q^^^'- are de Witt metrics defined as 

with inverse 

Qijki = -^{gikgji + giWjk) - -gijgu , Qijki = -{fikfji + fufjk) - -fijfki (2.9) 

that obey the relation 

Q^JklQ''''^'' = li^r^j + s?sT) ' ^u-fc^^'""" = li^r^j + s?sT) ■ (2-10) 

Using ( p.6| ) we rewrite the action ( p.l| ) into the form that is suitable for the Hamiltonian 
analysis 

S = j dtL = Ml j d^xdt^N[Kijg'^^^Kki + R^'^^] + 

+ M] j d^^dty/]M[kijg''^''^kki + R^f^] -II j d^xdtg^^'^f^/^VNMV . 

(2.11) 



■^We ignore the boundary terms. 
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Then from ( |2.11| ) we determine following conjugate momenta 



bL dL 



6dtN ' 6dtM 



(2.12) 



and then using the standard procedure we derive following Hamiltonian 



(2.13) 



where 



7^;^) = -25i,Vfc7r*=^' , 7^S^) = -2/,,Vfcp'=^' . 

(2.14) 

An important point is to identify four constraints that are generators of the diagonal 
diffeomorphism. In order to do this we proceed as in |l^ and introduce following variables 



N = VNM , n = \—, N' = -{N' + V), n 



M ' 2' ' ' ViVM ' 

N -. I . - . -. I . - 

N = Nn , M = — , V = N' - -n'N , N' = N' + -n'N , 
n ' 2 ' 2 ' 

(2.15) 

where again clearly their conjugate momenta are the primary constraints of the theory 

Pn-0 , Pn-0, , P^-0 . (2.16) 

Note that the canonical variables have following non-zero Poisson brackets 
{5.,(x),^'='(y)} = i(<5f5j + 5^5,')5(x-y) , {/,,(x), /'(y)} = i(5f ^ + 5^5^)5(x - y) , 

{iV(x),P^(y)} =5(x-y) , {n(x), P„(y)} = 5(x - y) , 
{iVXx),P,(y)} = S]6{^-y) , {n^(x),p,(y)} = 5i,5(x - y) . 



(2.17) 



Further, using ( p. 15 ) we find following form of the matrix i^t 

iV* 1 
= a + -^Vi , H% = -^Vi , 



H\ = -N'a - ^VkN^ + a\N^ + Nw' , 
N 



(2.18) 
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where 



Vi 



(2.19) 

The crucial point for the existence of four constraints that generate the diagonal diffeo- 
morphism is to show that the potential V does not depend on N and TV*. To do this we 
will argue that the matrix iJt is similar to the matrix A'^u defined as 



(2.20) 



Our arguments are as follows. It is known that matrices are characterized by their charac- 
teristic polynomials. The characteristic polynomial of an n x n matrix X is given by 

p{X) = det(X - AI) = (-1)"[A" + ciA"-i + caA"-^ + . . . + c„_iA + c„] , (2.21) 

where ci , C2 , . . . , c„ are expressed using the powers of the traces of X. Explicitly, if we 
introduce the notation 

tk = Tr(X'=) (2.22) 

we find that in case of 4 x 4 matrix the characteristic polynomial is determined by following 
coefficients 

ci = -ti , 

C2 = -t2) , 

1 o 1 1 

C3 = --4 + 2*1^2 - , 

1 4 1 2 1 1 2 1 



(2.23) 



Let us now determine tk for the matrix H% 

t, = TtH = H^'^ = a + a\ = A^^ , 

ts = TtH'-^ = Mf'^H^pHf^ = a^ + Sv'wia + 3via'jW^ + a' ja\a\ 
— A'^ A'' 



= a'^ + Aa^ViW^ + 2(i;ju;')^ + AaVio'jW^ + Avia' ^a\w^ + ja\a^ia\ 



\^^ \p A"" 



(2.24) 
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We see that the characteristic polynomials of -fft and A^,/ are the same which also implies 
that corresponding minimal polynomials are the same. Further, by presumption H'^y is 
diagonalizable which means that the minimal polynomial is a product of distinct linear 
factors. However since the minimal polynomial of H^y is the same as the minimal polyno- 
mial of A^u we immediately find that A.^i, is diagonalizable to the same diagonal matrix 
which also implies that and A^u are similar. In other words there exist the matrix Tt 
so that 

i/t = nA".(r"'n • (2-25) 

Then we find 

V{H) = V(A) (2.26) 
since by definition the potential is given as the traces of various powers of H^y. Now the 



equation (2.26) is the desired result which shows that V does not depend on N and and 
hence they appear in the action linearly. Using these calculations we find the Hamiltonian 
in the form 

H = j d^x{NlZ + N^iZi) , (2.27) 

where 

n = nni'^ + -ni^^ + -n'nl'^ - -n^n^^ + 

^ n ^ 2 * 2 * 

(2.28) 

Now we proceed to the analysis of the requirement of the preservation of the primary 
constraints ( 2.16| ) 



^tP^ = {Pi,H} = -TZ^^O , 

dtPn = {Pn, H] = -7^(^) + \n\!~^ - f^'g'/'f/^^^ ^Gn^O 

on 

dm = te,^} = + \j^P - ^^'9'"f'"^Z " " ° 



(2.29) 



As a result we have following total Hamiltonian 

Ht = j d^^{Nn + N'1li + VnPn + V'Pi + VnPn + v'pi + U^Gn + U^Gi) , 

(2.30) 

where V/v, V^,Vn,v'^,u'^, are Lagrange multipliers corresponding to the constraints Pfq « 
0, Pj « 0,pri ~ 0,pi 0, 0,^j Rs 0. In the next section we determine the algebra of 
constraints TZ^TZi. 
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3. Algebra of Constraints TZ, TZi 



For the consistency of the theory it is important to show that the constraints TZ and TZi 
are the first class constraints. To proceed it is useful to introduce the smeared form of the 
constraint IZ 

Tt{N) = [ d^xN{x)1Z{x) . (3.1) 



Instead of the constraint TZi we introduce the constraint TZi that is given as an extension 
of the constraints iZi by appropriate combinations of the primary constraints Pn,Pi 



TZi = TZ^f^ + TZ\f^ + diupn + diu^pj + dj{n^pi 
with equivalent smeared form 



(3.2) 



(3.3) 



Then using (|2.17D we determine the Poisson brackets between T5(A^*) and the canonical 
variables 



{TsiN'),gij} = 


-N^dkQij - diN^gkj - 




{T5(A^^),vr'^'} = 






{Ts{N%n,} = 


-N^dkfij - d.N^fkj - 




{TsiN^),p''} = 






{Ts{N'),n} = 


-N'diU , 




{Ts{N'),7Tn} = 


-^^{N'^ln) , 




{Ts(iV*),n*} = 


-N^dku' + djN'n^ , 




{T5(iV*),vr,} = 


-dkiN'm) - m^TTk . 





(3.4) 

With the help of these results we easily find how various components A^iy transform under 
spatial diffeomorphism 

{Ts{N'),a} = -N'd.a , 
{Ts{N'),Vi} = -N^djVi-m'vk , 



and also 

{t:s{N'),Ts{M^)] = TsiiN^djM' - M^djN')) . 
This result suggests that Ts{^^) is the generator of the spatial diffeomorphism. 



(3.5) 
(3.6) 
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For further purposes it is convenient to introduce the smeared forms of the constraints 



n\;^'^'^ and 7^K)'(^) 



Tf,(iV) = j c^3xiV(x)7^[,^^(x) , T^(iV) = j S^N{^)n\l\^) , 



d^xiV*(x)7lJ^^(x) , TUN' 



(3.7) 



It is well known that these smeared constraints have following non-zero Poisson brackets ^ 
{T?.(iV),Tf.(M)} = T^iNdiM - Mm)g'^) , 



{T^(iV),T^(M)} 
{T|(iV*),Tf.(M)} 
{T^(iV*),T^(M)} 
{T|(iV^),T|(M^)} 
{T^(iV^),T^(M^)} 



T^Ui^diM - Mm)P) , 
T^N'diM) , 
T^riN'diM) , 

TUiN^djM' - M^djN')) , 
T{{{N^djM' - M^djN')) . 



(3.8) 



As the next step we determine the Poisson bracket between Ts{N^) and Tj'(A^). We firstly 
determine following Poisson bracket 



1/4 fl/i 



6V 

5n 



kl 



6g 



(3.9) 



On the other hand we know that V depends on A'V through the trace. As a result we find 
that V depends on the eigenvalues Aj that are solutions of the characteristic polynomials. 
Further, we know that the coefficients of this polynomial are functions of ti, . . . ,^4. As a 
result if we show that ti transform as scalars under spatial diffeomorphism we find that Aj 
and consequently V transform as the scalar as well. In fact, using the explicit form of ti 
given m (|^) and also using the Poisson brackets ( |3.5D we easily find that ti transform as 
scalars under spatial diffeomorphism. These results imply following Poisson bracket 



{Ts{N),g'/'f/''v} = -N'dk[g'/'f/^V] - duN^g^l' f^/'V . 

Finally using the Poisson brackets ( |3.8| ) we find 

{Ts{N'),Tt{M)] = TriN'diM) , 
{Ts(iV*),G„(M)} = GniN^m) . 



(3.10) 



(3.11) 



'See, for example 
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In case of Qi the situation is more complicated since we have to determine the variation of 
ti with respect to and then its Poisson bracket with T5'(A^*). For example, in case ti 
and t2 we obtain 



6ti 

6t2_ 

dm} 



In the same way we can calculate 4% and -Pr and then we find 



5t^ 



on'' 



(3.12) 



(3.13) 



With the help of these results we easily obtain 



5V_ 

5n^ 



SV_ 

5n^ 



Finally using (gj) and ( pl4D we find 



(3.14) 



(3.15) 



Collecting all these results we find that TZi are the first class constraints whose smeared 
forms correspond to the generator of the diagonal spatial diffeomorphism. 

On the other hand more interesting is to determine the Poisson bracket between 



smeared forms of the Hamiltonian constrains (3.1). Following the same calculations as 
in we find 

{Tt{N),Tt{M)} = TsiiNdiM - MdiN)n^g'^) + Ts{{Nd^M - MdiN)\r) - 

- Gs{{NdiM - MdiN)n^g'^) - GriiNdiM - Md,N)n') + 

1 



+ GsiiNdiM - MdiN)^P) + / d^x{NdiM - MdiN)^'[V] , 



where we defined the smeared forms of the constraints Qi and Qn 

Gt{N) = [ d3xAr(x)g„(x) , Gs{N') = [ d^^N\x)gi{x) , 



(3.16) 



(3.17) 



and where S*[V] is defined as 



2 ^i 6V 1 



6V 



-n'g'^ 



—n - -n n- 



(3.18) 
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We see that the Poisson bracket between the smeared forms of the Hamiltonian constraints 
( |3.16[ ) vanish on the constraint surface on condition that S*[V] is zero. In order to calculate 
S*[V] we proceed in the similar way as in case when we calculated the Poisson bracket 
between Tsi^^) ^-nd V. Explicitly, we know that A^iy is similar to diagonal matrix with 
eigenvalues Ai, A2, A3, A4. In other words we have 



V(A) = V(A,) . (3.19) 

On the other hand we know that Aj are solutions of the characteristic polynomial when the 
coefficients are functions of ti , . . . , t4 so that 

Xi = X,{h,...,U) . (3.20) 

=0 ,i = 1,2,3,4 (3.21) 
S*[Aj] = (3.22) 



Then if we show that 



we find that 



and consequently 



S^[V(A)] = S^[V(A)] =0 . (3.23) 



We have already found in [^] that = ^^[^2] = 0. The case of and is more 
intricate. For example, the explicit form of ts is equal to 

- ^n'f.kg^'fimn^ + 4(n7.fc/7^n^n™)(n7rsn^) - 

(3.24) 

Then using the explicit form of and after some calculations we find 

E'[t3]=0. (3.25) 
In the same way we find S*[t4] = 0. In summary we have 

= , for j = 1,2,3,4 . (3.26) 



With the help of this result we obtain the fundamental result that the Poisson bracket 
between Hamiltonian constraint ( [j.l6 ) vanishes on the constraint surface. This result 
together with ( |3.6| ) and ( 3.16| ) shows that Tr(A^) and T5(A^*) are the first class constraints 
that are generators of the diagonal diffeomorphism. 
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4. Time Evolution of Remaining Constraints 



In this section we analyze the time evolution of the constraints PmPi,Gn,Gi- Note that the 
total Hamiltonian is given in (Ul). For further purposes we calculate following Poisson 
brackets 



{Pn(x),^i(y)} = fl^ 



(x)(5(x - y) = A„j(x)J(x - y) , 



{pn(x),^„(y)| = -^t/„,(x)(5(x - y) + 
+ n(x) 

n(x) 



5^V 



6n^6n 
(x 



n 

1 ^ n 

n n 

sy_ 

5n 



2„l/4fl/4^^ ^2^1/4^1/4^ 



n 



5n 

1/4^1/4^ _ 

52 



(5n(x)(5n(y) 



5(x-y) 



52 



n 



^2^1/4/1/4 



<5(x-y) 
(52V 



n*— ^ - V + n— + n^^TT- 
on' on o'^n 



(5(x - y) = A„„(x - y) 



te(x),an(y)} = ^2gi/4/i/4^(x)^-(x _ y) ^ A,„(x)5(x - y) 

drv-on 



(4.1) 



We introduce following common notation Pa = (PmPi) and Ga = {GniGi)- Then the 
requirement of the preservation of the primary constraints pa takes the form 



dtPa = {Pa^Hr} = Gn + j d^^U^ {^) {Pa, Gb{^)} 
(i3xn^(x) {pa,Gbi^)} = Aabvf' = . 



(4.2) 

We have four equations for four unknown u"". We have to distinguish two cases. The case 
when det 7^ and the case when det /\ab = 0- 

4.1 The case det A^;, / 

This situation is particular simple and we believe that this is in fact the situation in all 
examples of the bimetric theories of gravity that are invariant under diagonal diffeomor- 
phism. In this case we find that the solution of the equations (4.2) is v!' = 0. Then we can 
proceed to the analysis of the preservation of the constraints Ga 



dtGa = {Ga, Ht} ~ J d'j^ (iV(x) {Ga, i?(x)} + v' {Ga,Pbi^)}] 

= J d^xNix) {Ga,R{^)} - AabV^ = 



(4.3) 
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that again using the fact that det A^fe 7^ can be solved for f " as functions of the canonical 
variables and A^. In other words we find that Qa,Pa are the second class constraints that can 
be solved for n, n* at least in principle. Then four first class constraints Pf^,Pi can be gauge 
fixed and hence we can eliminate Pjy, Pj and N,N^ as dynamical variables. Then remaining 
four first class constraints 1Z,TZi can be again gauge fixed so that we can eliminate eight 
degrees of freedom from 24 degrees of freedom gij , tt^-^ , fij , p^-^ so that we have 16 physical 
degrees of freedom where 4 correspond to the massless graviton, 10 to the massive one 
and 2 corresponding to the scalar mode at least at the linearized level. Of course, this 
Hamiltonian analysis cannot answer an important question which is the identification of 
the physical metric that couples to the matter fields, for recent detailed analysis, see |^. 

To conclude this section we have to mention the form of the symplectic structure 
on the reduced phase space spanned by = (gij , vr*-' , fij , p*-' ) . Due to the fact that we 
have eight second class constraints Pn,Pi,Gn, Qi we should replace the Poisson brackets 
with corresponding Dirac brackets. Explicitly, let us introduce the common notation for 
all second class constraints XA = iPa,Ga)- Then we find following matrix of the Poisson 
brackets of the second class constraints XA 

r / N / ^ f A(x)(^(x-y)\ 

{XA{^),XA{y)}^eAB= [ ..... . ^'^ /M , 4.4 

\-A(y)5(x-y) J^(x,y) J 

where the matrix il(jb(x,y) is defined as 

{ga{^),Gb{y)} = ^ab{^,y) ■ (4.5) 

Since A is invertible matrix we can find the matrix inverse to Qab in the form 

(a-i)AB^ ( A-i(x)l^(x,y)A-i(y) - A-i(x)5(x - y) 
^ ^ A-i(y)5(x-y) 

Following the standard analysis of the system with the second class constraints we replace 
the Poisson bracket by corresponding Dirac bracket defined as 

{F,G}^ = {F,G}- J d3xd3y{i^,XA(x)}(G-i)^^(x,y){xB(y),G} . 

(4.7) 

Now we see that for the variables on the reduced phase space spanned by = {gij , vr*-' , fij , p*-^ ) 
the Dirac bracket coincides with the Poisson bracket. Explicitly, we have 

{z"(x),z/^(y)}^ = {z"(x),z^(y)}- 

- I d3zdV{^"(x),XA(z)}(e-i)^^(z,z'){xB(z'),^''(x)} = {z"(x),z/^(y)} - 

- / dVz'^ ^{z"(x),XA(z)}(e-i)^^(z,z'){xB(z'),^''(x)} = {z"(x),z/^(y)} 

(4.8) 

due to the fact that (B-^)^-^ = for ^, 5 = 3, . . . , 7 are zero. 
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4.2 The case det Aab = 

In this case it is possible to find the vector such that 

^ab4 = . (4.9) 
Then we define new constraints P, Q defined as 

p = ngp„ ,a = ^^Sg„ (4.10) 

together with 

Pa = Pa — ^^OaP , Qa = Qa ^—^UQaG , Uoa = SabUo (4.11) 

that by definition obey the conditions 

u^Pa = 0, nga„ = (4.12) 

which means that we have 6 independent constraints. As a result we consider the total 
Hamiltonian in the form 

Ht = J d^xiNTZ + N'iii + V^vPjv + V'Pi + vP + v'^Pa + uQ + U^Ga) , 

(4.13) 

where 

u"'UOa ^io^Oa 

where v°',u°' are arbitrary four dimensional vectors and we clearly see that v"',u°' are 
orthogonal to Uq so that they have 3 independent components. 

Now we proceed to the analysis of the time evolution of the constraints P, Pa,G,Ga- 
Note that the constraint P has vanishing Poisson bracket with G and Ga since 

{P(x),g~(y)} = u^A,,45{^-y) = , 
{P(x), Ga{y)} = 4Ab„5(x - y) = 

(4.15) 

so that it is trivial to sec that P is preserved during the time evolution of the system and 
that it is the first class constraint. Let us now analyze the time evolution of the constraint 



Pa 



using 



dPa = {Pa, Ht] ~ J d^^{u{^) [Pa, G{^)} + u\^) {p„, Gb{^)]) « A^bU^ = 

(4.16) 

{P„(x),g(y)} = {Pa{^),Gb{y)}4 = ^ab4S{^ - y) = . (4.17) 
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Now due to the fact that u'^ is orthogonal to Uq we find that the only solution of ( 4. IE ) is 
given by = 0. 

We proceed to the analysis of the time evolution of the constraint G,Qa- We start with 
the constraints Ga 

dtga = {ga,HT}- J d'x (Af(x){g,,^(x)}+n(x){g,,g(x)}) +^^Ab„ = 

(4.18) 

that again using the definition of given ( 4.14| ) can be solved for v"" as function of the 
canonical variables and and u. This result implies that Pa,Qb are the second class 
constraints. 

Now we come the most intricate analysis which is the time preservation of the constraint 

Q 

dtg = [g,HT] ~ j d3x(iv(x){g,^(x)} +n(x){g,^(x)}) = o . 

(4.19) 

Note that generally |^(x),^(y)| ^ 0. However we see that we have to analyze the time 
evolution of the constraint TZ as well ^ so that 

dtU = {n, Ht] ~ j (i^xn(x) j:^, g(x)} = . 

(4.20) 

To proceed further we have to determine the Poisson bracket between TZ and Q. This is very 
complicated due to the complex form of these constraints. As an example we determine 
the Poisson bracket between T^t{N) and Gt{M) 

{Tt{N),Gt{M)} = -T^g{n{NdiM - MdiN)g'^ - MNdiug'^) + 

+ T^(n"^(A^9iM - MdiN)p - n'^dinNMp) - 

- ]^T9^{n'NdiM) - ]^T{,{n'Ndi{n-^M)) + 

+ 2l? [ d^xMN I — mr'^^gnkiG''^ J-^p'^^gkiaF'^ ] + 

^ J \M^V9 ' Mfn^'^ ^'"''^ J ^ 

(4.21) 



*Note that this was not necessary in case when det Aat since in this case we found that u"' — and 
hence the constraint TZ is trivially preserved. 
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where 



^9kl 



Sfki 



kl 



(4.22) 



In the same way we calculate the Poisson bracket between T(A^) and Gs{N^) and we derive 
similar result. Then in principle we could express the Poisson bracket between T(A^) and 
J (i^xM(x)^(x) in the form 

Tt{N),I d3xM(x)g(x)| = 

d^z{N{z)M{z)F{z) + ^,^N{z)Y'{z)M{z) + N{z)^,^M{z)W'{z)) , 



(4.23) 

where the explicit form of F, V*, W* can be derived from the Poisson brackets between 
T(A^) and Gt{N),Gs{N^) and from the known vector ng. For our purposes it is useful 
the local form of the Poisson bracket ( 4.23| ) 



d 



{^(x), g(y)} = 5(x - y)F(x) + ^6{x - y) V(y) + ^5(x - y)W^(x) . (4.24) 



With the help of this result we find that ( 4.201) takes the form 



dtU = n(F - diW') + -^(W* - V* 
ax* 



0. 



(4.25) 



In order to find the meaning of this result we can argue in the same way as in |28[| where 
more details can be found. Briefly, as we can deduce from ( 4.21| ) we have that W 7^ V* so 



that we should interpret the equation ( 4.25 ) as the differential equation for it(x) and this 
equation could be solved for u(x) at least in principle. Further, using ( [4.24| ) we find that 



(4.19) becomes differential equation for A^(x) 



dN 



dtQn = iV(-F + diW) + —^{W 



V) + 



+ 



j c^V^i(y){an(x),an(y)} = . 



(4.26) 



This equation can be now solved for A^(x) as function of the canonical variables since the 



Lagrange multiplier u has been already determined by the equation ( 4.20 ). These results 
imply that 1Z and Q are the second class constraints. On the other hand we could say 
that by consistency the theory should have the first class Hamiltonian constraint as a 
consequence of the fact that the theory is invariant under diagonal diffeomorphism. In 
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other words we are tempting to say that in bimetric theories that are invariant under 
diagonal diffeomorphism the case when det Aab = cannot occur ^ . 



4.3 Square root potential 

Recently it was shown that the non-linear massive gravity with the square root form of 
the potential is free from ghosts. The extension to the bimetric gravity was performed 
in ||2^ where it was also argued that resulting theory should be ghost free. This fact was 
questioned in |28| where the Hamiltonian analysis of the bimetric gravity introduced in |25| 



was performed. We argued that the theory with redefined shift function lost the manifest 
diffeomorphism invariance and so that it difficult to find the generator of the diagonal 
diffeomorphism. We also argued that the number of the physical degrees of freedom is 16 
so that the scalar mode cannot be eliminated. 

In order to resolve the apparent issue of the impossibility to identify the generators 
of the diagonal diffeomorphism in the non-linear bimetric theory of gravity with redefined 
shift function we can certainly use the analysis of the general bimetric theory of grav- 
ity presented in this paper. Obviously we can identify four first class constraints IZ^lZi 
corresponding to the diagonal diffeomorphism. On the other hand we can ask the question 
whether there is possibility to identify an additional constraint that could eliminate the 
scalar mode. To be more concrete, let us consider following square root potential 



V = Try^ = TV\/a^ • (4.27) 
The problem is that it is not easy task to find the square root of the matrix A^^, since we 



were not able to find such redefinition of the canonical fields as in 14 1 that could allow 
us to find the explicit form of the square root of A'^u- For that reason we try to find such a 
matrix in the following way. We temporally write with the parameter e so that n* — )• en* 
when we take e = 1 in the end. Then we can write 



where 



A = eOE(o) + eE(i) + e^^^"^^ + e^E^^) , (4.28) 



\ / ^ 







(4.29) 



^The situation is different in the miracle case when W = V . In this case it is more natural to interpret 
g(3) = F _ 9iV' = as the new constraint rather than impose the condition ii = Q [Q. Then of course 
the equation (4.19) vanishes on the constraint surface C?'^' ~ and the requirement of the preservation 



of the constraint t/'"^' ~ leads to the emergence of the new constraint ~ on condition that the 
Poisson bracket between t/''^' and TZ does not contain the derivative of the delta function. However due to 
the complex form of the constraints TZ and Qa and corresponding Poisson brackets we believe that such a 
case is highly improbable. 
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We can solve v A as the expansion with respect to e 

oo 



(4.30) 



n=0 



Taking square of given expression and comparing expressions of given order in e we find 
E(o) = b^bW , 

+ b(2)b(0)+B«B« , 
+ B«b(2)+b(2)b« + b(3)bW , 

+ b(2)b(2) + b(3)b(i) + bWb(4) + b(4)b(o) , 



e(i) = b(0)b(i 

e(2) = b(0)b(2 

e{3) = b(°)b(3 

= b(i)b(3 



(4.31) 

From the first equation in ( [4.31| ) we find 

bW = ( nf M , (4.32) 




where 



7\ = \/9''fkj, 7\l'j=9"'fkj . (4.33) 



Further, the second equation in ( 4.31D can be solved as 

B'" = (,„_,„ , ''] (4.34) 




where 

M\ = ^S\ + f, . (4.35) 

Repeating this procedure we could determine all matrices B^") at least in principle. On 
the other hand we see that the resulting expression contains all powers of n and n*. Then 
from (4.1) we can deduce that Aab 7^ and hence Pa,Ga are the second class constraints 
without possibility to eliminate additional scalar mode. 



5. Conclusion 

Let us outline the results derived in this paper. We performed the Hamiltonian analysis 
of bimetric theory of gravity with general form of the potential. We found the first class 
constraints that are generators of the diagonal diffeomorphism. Then we analyzed the 
requirement of the preservation of remaining constraints during the time development of 
the system. We show that there are two possibilities. In the first case when the determinant 
det Aab 7^ we found that there are eight second class constraints. Unfortunately these 
second class constraints are not sufficient for the elimination of the scalar mode. On the 
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other hand we showed that in the second case det Aab = we find the first class constraint 
P. However then we argued that in this case the Hamiltonian constraint TZ together with 
the constraint Q are the second class constraints. On the other hand since we know that 
the theory is invariant under diagonal difi'eomorphism we suggests that the Hamiltonian 
constraint should be the first class constraint and consequently the case when det Aab = 
should not occur. All these results could also imply that the ghost mode is still presented 
in the non-linear bimetric theory of gravity. 

As the extension of this work we suggest the Hamiltonian analysis of the multimetric 



theories of gravity that were introduced in paper ||31[] and we hope to perform such analysis 
soon. 

Acknowledgement: This work was supported by the Grant agency of the Czech republic 
under the grant P201/12/G028. 

References 

C. J. Isham, A. Salam and J. A. Strathdee, "F-dominance of gravity," Phys. Rev. D 3 (1971) 
867. 

S. F. Hassan, A. Schmidt-May and M. von Strauss, "Bimetric Theory and Partial 
Masslessness with Lanczos- Lovelock Terms in Arbitrary Dimensions," arXiv:1212.4525 
[hep-th]. 

S. F. Hassan, A. Schmidt-May and M. von Strauss, "On Consistent Theories of Massive 
Spin-2 Fields Coupled to Cravity," arXiv:1208.1515 [hep-th]. 

K. Nomura and J. Soda, "When is Multimetric Cravity Chost-free? ," Phys. Rev. D 86 
(2012) 084052 [arXiv:1207.3637 [hep-th]]. 

V. Baccetti, P. Martin-Moruno and M. Visser, "Null Energy Condition violations in bimetric 
gravity," JHEP 1208 (2012) 148 [arXiv:1206.3814 [gr-qc]]. 

V. Baccetti, P. Martin-Moruno and M. Visser, "Massive gravity from bimetric gravity, " 
Class. Quant. Grav. 30 (2013) 015004 [arXiv:1205.2158 [gr-qc]]. 

M. Banados and S. Theisen, "Three-dimensional massive gravity and the bigravity black 
hole," JHEP 0911 (2009) 033 [arXiv:0909.1163 [hep-th]]. 

M. Banados, A. GomberofF, D. C. Rodrigues and C. Skordis, "A Note on bigravity and dark 
matter," Phys. Rev. D 79 (2009) 063515 [arXiv:0811.1270 [gr-qc]]. 

D. Bias, C. DefFayet and J. Garriga, "Clobal structure of bigravity solutions," Class. Quant. 
Grav. 23 (2006) 1697 [hep-th/0508163]. 



[10 



[11 



[12' 



[13 



T. Damour and I. I. Kogan, "Effective Lagrangians and universality classes of nonlinear 
bigravity," Phys. Rev. D 66 (2002) 104024 [hep-th/0206042]. 

C. de Rham and G. Gabadadze, "Generalization of the Fierz-Pauli Action," Phys. Rev. D 82 

(2010) 044020 [arXiv: 1007.0443 [hep-th]]. 

C. de Rham, G. Gabadadze and A. J. Tolley, "Resummation of Massive Cravity," Phys. Rev. 
Lett. 106 (2011) 231101 [arXiv:1011.1232 [hep-th]]. 

S. F. Hassan and R. A. Rosen, "On Non-Linear Actions for Massive Gravity," JHEP 1107 

(2011) 009 [arXiv:1103.6055 [hep-th]]. 



-18- 



[14] S. F. Hassan and R. A. Rosen. "Resolving the Ghost Problem in non-Linear Massive 
Gravity," Phys. Rev. Lett. 108 (2012) 041101 [arXiv: 1106.3344 [hcp-th]]. 

[15] S. F. Hassan, R. A. Rosen and A. Schmidt-May, "Ghost-free Massive Gravity with a General 
Reference Metric," JHEP 1202 (2012) 026 [arXiv:1109.3230 [hep-th]]. 

[16] J. Kluson, "Note About Hamiltonian Structure of Non-Linear Massive Gravity," JHEP 1201 
(2012) 013 [arXiv:1109.3052 [hep-th]]. 

[17] J. Kluson, "Hamiltonian Analysis of 1+1 dimensional Massive Gravity," Phys. Rev. D 85 
(2012) 044010 [arXiv:1110.6158 [hep-th]]. 

[18] A. Golovnev, "On the Hamiltonian analysis of non-linear massive gravity, " Phys. Lett. B 
707 (2012) 404 [arXiv:1112.2134 [gr-qc]]. 

[19] J. Kluson, "Comments About Hamiltonian Formulation of Non-Linear Massive Gravity with 
Stuckelberg Fields," JHEP 1206 (2012)170 [arXiv:1112.5267 [hep-th]]. 

[20] J. Kluson, "Remark About Hamiltonian Formulation of Non-Linear Massive Gravity in 
Stuckelberg Formalism," Phys. Rev. D 86 (2012) 124005 [arXiv:1202.5899 [hep-th]]. 

[21] J. Kluson, "Non-Linear Massive Gravity with Additional Primary Constraint and Absence of 
Ghosts," Phys. Rev. D 86 (2012) 044024 [arXiv: 1204.2957 [hep-th]]. 

[22] J. Kluson, "Note About Hamiltonian Formalism for General Non-Linear Massive Gravity 
Action in Stuckelberg Formalism," arXiv: 1209.3612 [hep-th]. 

[23] S. F. Hassan and R. A. Rosen, "Confirmation of the Secondary Constraint and Absence of 
Ghost in Massive Gravity and Bimetric Gravity," JHEP 1204 (2012) 123 [arXiv:1111.2070 
[hep-th]]. 

[24] S. F. Hassan, A. Schmidt-May and M. von Strauss, "Proof of Consistency of Nonlinear 
Massive Gravity in the Stuckelberg Formulation, " Phys. Lett. B 715 (2012) 335 
[arXiv:1203.5283 [hep-th]]. 

[25] S. F. Hassan and R. A. Rosen, "Bimetric Gravity from Ghost-free Massive Gravity," JHEP 
1202 (2012) 126 [arXiv:1109.3515 [hep-th]]. 

[26] E. Gourgoulhon, "3-1-1 formalism and bases of numerical relativity," gr-qc/0703035 [GR-QC]. 

[27] R. L. Arnowitt, S. Deser, C. W. Misner, "The Dynamics of general relativity," 
[gr-qc/0405109]. 

[28] J. Kluson, "Ls Bimetric Gravity Really Ghost Free?," arXiv: 1301. 3296 [hop-th]. 

[29] S. A. Hojman, K. Kuchar and C. Teitelboim, "Geometrodynamics Regained, " Annals Phys. 
96 (1976) 88. 

[30] J. Kluson, "Hamiltonian Formalism of Particular Bimetric Gravity Model," arXiv: 121 1.6267 
[hep-th]. 

[31] K. Hinterbichler and R. A. Rosen, "Interacting Spin-2 Fields, " JHEP 1207 (2012) 047 
[arXiv: 1203.5783 [hep-th]]. 

[32] V. O. Soloviev and M. V. Tchichikina, "Bigravity in Kuchar's Hamiltonian formalism. 2. 
The special case," arXiv:1302.5096 [hep-th]. 



- 19 - 



